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1. Introduction 



Stochastic partial differential equations (SPDEs) combine the features of partial differential equa- 
tions and Ito equations. Such equations play important roles in many applied fields such as the 
filtering of partially observable diffusion processes, genetic population and other areas. We study 
the following stochastic partial differential equation (in short SPDE) for a system-valued of pre- 
dictable random field Ut (x) = u {t, x) , satisfying the following equation: 



dut{x) + [Cut{x) +f{t, X, ut{x),'\/ut(j{x))) dt + g{t, x, ut{x),Vuta{x)) ■ dBt = 0, (1-1) 

over the time interval [0,T], with a given final condition ~ ^ and non-linear deterministic 
coefficients /, g. Lu = {Lui, • • • , Luk) with L a second order differential operator. The differential 
term with dBt refers to the backward stochastic integral with respect to a Z-dimensional Brownian 
motion on J^, P, {Bt)t>o) ■ We use the backward stochastic integral in the SPDE because we will 
employ the framework of Backward Doubly Stochastic Differential Equation (BDSDE) introduced 
first by Pardoux and Peng [26]. They have given a probabilistic representation for the classical 
solution ut{x) of the SPDE (1.1) (written in the integral form) in term of the following class of 
BDSDE's: 

y,*'" = $(X^") + y f{r,X','''X'":K'")dr + J g(r,X;'^F/•^Z*'")dB,-y (1.2) 

where {Xl'^)t<s<T is a diffusion process starting from x at time t driven by the finite dimensional 
brownian motion (Wt)t>o a-nd with infinitesimal generator L. More precisely, under some regularity 
assumptions on the final condition $ and coefficients / and g , they have proved that ut{x) — Y^'^ 
and Vut(T{x) = Z^"" , V(t,x) G [0,r] x . Then, Bally and Matoussi [5] (see also [22] ) showed 
that the same representation remains true in the case when the final condition (respectively the 
coefficients / and g) is only measurable in x (resp. are jointly measurable in (t, a;) and Lispchitz 
in u and Vit). In this paper, weak Sobolev solution of the equation (1.1) has been considered, and 
the approach was based on stochastic fiow technics (see also [17, 18]). Moreover their results were 
generalized in [22] in the case of a larger class of SPDE's (1.1) driven by a Kunita-Ito non-linear 
noise (see [17, 18, 19] for more details). In particular, the Kunita-Ito non- linear noise covers a class 
of infinite dimensional space-time colored-white noise (see [12], [27], [15]). Generally, the explicit 
resolution of semi-linear SPDEs is not possible, so it is then necessary to resort to numerical 
methods. 

The first approach used to solve numerically nonlinear SPDEs is analytic methods, based on time- 
space discretization of the SPDEs. The discretization on space can be achieved either by finite 
differences, or finite elements and spectral Galerkin methods. But most numerical works on SPDEs 
have concentrated on the Euler finite-difference scheme. Gyongi and Nualart [13] have proved that 
these schemes converge, and Gyongy [14] determined the order of convergence. J.B. Walsh [28] 
investigated schemes based on the finite elements methods. He studied the rate of convergence of 
these schemes for parabolic SPDEs, including the Forward and Backward Euler and the Crank- 
Nicholson schemes. He found substantially similar rate of convergence to those found for finite 
difference schemes. The spectral Galerkin approximation was used by Jcntzcn and Kloedcn [15]. 
They based their method on Taylor expansions derived for the solution of the SPDE, under some 
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regularity conditions. 

Lototsky, Mikulevicius and Rozovskii in 1997 [20] used the spectral approach for the numerical 
estimation of the conditional distribution solution of a linear SPDE known as the Zakai equation. 
Further developments on spectral methods can be found in Lototsky [21]. Very interesting results 
have been obtained by Gyongy and Krylov [12] where they considered a symmetric finite difference 
scheme for a class of linear SPDE driven by infinite dimensional Brownian motion. They have 
proved that the approximation error is proportional to fc^ where k is the discretization step in 
space and by the Richardson acceleration method they have even got the error proportional to k^. 

The other alternative for resolving numerically SPDEs is the probabilistic approach. It requires 
weaker assumptions on the SPDE's coefficients. In the deterministic PDE's case i.e. 5 = 0, the nu- 
merical approximation of the BSDE has already been studied in the literature by Bally [3] , Zhang 
[29], Bouchard and Touzi [6], Gobet, Lemor and Warin[ll] and Bouchard and Elie [7]. Zhang [29] 
proposed a discrete-time numerical approximation, by step processes, for a class of decoupled FBS- 
DEs with possible path-dependent terminal values. He proved an L^-type regularity of the BSDE's 
solution, the convergence of his scheme and he derived its rate of convergence. Bouchard and Touzi 
[6] suggested a similar numerical scheme for decoupled FBSDEs. The conditional expectations 
involved in their discretization scheme were computed by using the kernel regression estimation. 
Therefore, they used the Malliavin approach and the Monte carlo method for its computation. 
Crisan, Manolarakis and Touzi [8] proposed an improvement on the Malliavin weights. Gobet, 
Lemor and Warin in [11] proposed an explicit numerical scheme. In the case when g 7^ and it 
does not depend on the control variable z, Aman [1] proposed a numerical scheme following the 
idea used by Bouchard and Touzi [6] and obtained a convergence of order h of the square of the 
L^- error {h is the discretization step in time). Aboura [2] studied the same numerical scheme 
under the same kind of hypothesis, but following Gobet and al [10]. He obtained a convergence of 
order \/h in time and used the regression Monte Carlo method to implement his scheme, following 
always [10]. 

In our work, we extend the approach of Bouchard- Touzi-Zhang in the general case when g de- 
pends also on the control variable z. We wish to emphasize that this generalization is not obvious 
because of the strong impact of the backward stochastic integral term on the numerical approxi- 
mation scheme. It is known that in the associated Stochastic PDE's (1.1), the term g{u, Vu) leads 
to a second order perturbation type which explains the contraction condition assumed on g with 
respect to the variable z (see [26], [24]). Our scheme is explicit in Y and implicit in Z. We prove 
the convergence of our numerical scheme and we give the rate of convergence. The square of the 
L^- error has an upper bound of order the discretization step in time. As a consequence, we get 
a numerical scheme for the weak solution of the associated semi linear SPDE. We give also a rate 
of convergence result for the later weak solution. Then, we propose a numerical scheme based on 
iterative regression functions which are approximated by projections on vector space of functions 
with coefficients evaluated using Monte Carlo simulations. Finally, we present some numerical 
tests. Compared to the deterministic numerical method developed by Gyongy and Krylov [12], the 
probabilistic approach could tackle the semilinear SPDE which could be degenerate and needs less 
regularity conditions on the coefficients than the finite difference scheme. However, the rate of con- 
vergence obtained (as the classical Monte Carlo method) is clearly slower than the results obtained 
by difference and finite element schemes, but of course more available in higher dimension. 
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This paper is organized as follows. In section 2 we introduce preliminaries and assumptions. In 
Section 3, we describe the approximation scheme for the BDSDE and we state an upper bound 
result for the time discretization error. In section 4, we give a Malliavin regularity result for the 
solution of our Forward-Backward Doubly SDE's. Then, we show a L^-regularity result for the 
Z-componcnt of the solution of the BDSDE (1.2) which is crucial to obtain the convergence and 
the rate of convergence of our numerical scheme. Section 5 is devoted to the numerical scheme 
of the SPDE's weak solution. In section 6, we test statistically the convergence of this scheme by 
using a path dependent algorithm based on the regression Monte Carlo Method. Finally, we give 
some technical results in the Appendix. 

2. Preliminaries and notations 

2.1. Forward Backward Doubly Stochastic Differential Equation 

We consider a probability space supporting {IV^iO < s < T} and {-B^iO < s < T} two mutually 
independent standard Brownian motion processes, with values respectively in R'^ and in M' where 
T > is a fixed horizon time. 

We shall work on the product space fl := flw x where flw is the set of continuous functions 
from [0,r] into R'' and Hb is the set of continuous functions from [0,r] into R'. We fix t e [0, T]. 
For each s ^ [t,T], we define 

where T}^^ = -Wt,t < r < s}, and J'f = a{Br - Bs, s < r < T}. We take = T^j., 

= T^j, 'And T ^ T^' V . 
We define the probability measures Pw on {nw,J'^) and Pb on (17b, J"^). We then define the 
probability measure P := Pw 'E) Pb on (fi, x T^). Without loss of generality, we assume that 
J-^ and are complete. We denote by {^l,J^,P) our probability space. 

Note that the collection {J-l,s S [^,2^]} is neither increasing nor decreasing, and it docs not 
constitute a filtration. 

Given C > 0, we consider two functions b and a satisfying the Lipschitz condition i.e. for all 
X g R"*, x' e R"*, we have 

(HI) \b{x) ~ bix')\ + \\a{x) ~ a{x')\\ < C\x ~ x'\. 

We consider the following stochastic differential equation: 

dX*^^ = 6(X*'^)ds + cr{Xl^'=)dWs, s G [t, T], X*^^ = X, s< t. (2.1) 

We omit the dependance of the forward process X in the initial condition if it starts at time i = 0. 
We fix i < si < 52- For some real number p>2 and for any n e N, let Hfj([si, S2\) denotes the set 
of (classes of dP x dt a.e. equal) n dimensional progressively measurable processes u € [si, S2]} 
satisfying : 

W II^IIh-„([.i,.=]) ■■-=E[j:^\^u\Pdu]<oo, 
(ii) ipu is J^*-measurable, for a.e. u G [si,S2]- 

We denote similarly by §^([si, S2]) the set of continuous n dimensional processes satisfying : 
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(ii) ipu is J^^-measurable, for any u G [si, 52]- 

In the sequel, we shall omit the subscript [si, S2] in these notations when [si, S2] = [0, T]. We define 
the following functions : 



/ : [0, T] X X K'' X 
5: [0,T] X 



We shall make use of the following assumption : 

(H2) there exist two constants K > Q and < a < 1 such that for any 
(ti,a:i,yi,Zi),(t2,X2,2/2,^2) e [0,r] x M'' x R'^' x R^'x^ 

(i) |/(ti,a;i,yi,zi) - /(t2, 2:2, y2, ^2)! < K[y/\ti - + - a::2| + \yi - ^2! + - Z2\\), 

(ii) ||.g(ti,a;i,2/i, zi) - 3(^2,2:2,2/2,^2)11^ < -^'(Ki - ^2! + |a;i - 2:2 P + \yi - ^2!^) + c?\\zi - Z2II' 

(iii) |$(a;i)-$(a;2)| <is:|a;i-a;2|, 

(iv) sup 1/(^,0,0,0)1 + 0,0, 0)11 < if. 

0<t<T 

We consider the following backward doubly stochastic differential equation. For all t < s < T, 

= <I>(X^--). ^ 

We note that the integral with respect to (i3s,i < s < T) is a "backward Ito integral" and the 
integral with respect to {Ws,t < s < T) is a standard forward Ito integral. 

Pardoux and Peng [26] proved that there exists a unique solution (F, Z) € §fc([i, T]) x IHI^^^([i, T]) 
to the BDSDE (2.2). 




Remark 2.1 The contraction condition < a < 1 is needed to prove the existence and the 
uniqueness results for the BDSDE's solution (see Pardoux and Peng [26]) . 

From [9], [26] and [16], the standard estimates for the solution of the F-BDSDE (2.1)-(2.2) hold 
and we have the following theorem : 

Theorem 2.1 For any p > 2 there exists a positive constant C such that : 



E[ sup \Xr\^]<Cil + \x\'), 

t<s<T 



(2.3) 



E 



E 



sup 

Si < U<S2 



sup 

t<s<T 



|n*'n^+(y^' WzrfdsY^'] <C(l + |x|^), 
I^M _ yM|p] < c((l + )|52 - si|5 + ||Z*--||h^,,[.,...]), 



(2.4) 



(2.5) 



We denote by CI;{Rp,R''), C^{W,m) the set of functions of class C'' from Rp to R« whose partial 
derivatives of order less or equal to k are bounded and the set of functions of class C°° from R^' to 
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whose partial derivatives are bounded. For later use, we assume 

(H3(i)) $ e Cl{W^,ml')J £ Cl([Q,T] X M'' X M*^ X R'^^'^',R'=) 

and g £ Cl([Q,T] x K'^ x K'' x M'^^^ M'^''^') 

(H3(ii)) $ G CI{W^,'B!')J e Cfo([0,r] X R'^ X M'^ X R'^^'=,R''') 

and g G C^2([0,r] X R'' X M*^ X M'^^'^.M'^^') 

2.2. Numerical Scheme for decoupled Forward- BDSDE 

In order to approximate the solution of the BDSDE (2.2), we introduce the following discretized 
version. Let 

TT : to = < ti < ... < ijv = (2.6) 

be an equidistant partition of the time interval [0,r] i.e. h = ^ and i„ = nh, < n < N . 
Throughout the rest, we will use the notations AWn = Wt^^^ — Wt^ and AB„ = Bt^j_^ — Bt^, for 
n=l,...,N. 

The forward component X will be approximated by the classical Euler scheme : 



X^l = Xl_^ + b{Xl_^){U - t„-i) + a{Xl_J{Wt„ ~ J, for n = 1, ...,iV. 



(2.7) 



It is known that as N goes to infinity, one has sup E\Xt^-X^\'^ -> 0. 

0<n<A' 

Quite naturally, the solution {Y, Z) of (2.2) is approximated by {Y^ , Z^) defined by : 

Y,l = ^{X^), (2.8) 

and for < n < iV - 1, 

y,^ = Et^ [Yi^J + hEt^ [/(t„, C )] + [.g(t„+i, e,^+i)Ai3„], (2.9) 



hZ^" = Et. 



Y.'l^^AW: + .g(i„+i, e^+i)AS„AM^„* 



(2.10) 



where 

:= (x/:_„>;^z,^_J, e^:=(X,^,y,^£;,„[Z,^_J),Vn==l,..,iV. 
* denotes the transposition operator and Et^ denotes the conditional expectation over the cr-algebra 

We define also for aU n^O,..,N - 1, (F^, Z^)t^<s<t„+^ as the solution of the following BDSDE: 

I dY^ = -f{t^X)ds-g{tn+uQ^+^)tBs + ZfdWs. ^^.u) 
1 ^t^+i given by our numerical scheme. 

This is the continuous approximation of the solution of the BDSDE (2.2). 
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Remark 2.2 For the approximation of , (2.10) is well-defined, indeed if we set 



= 



i > 1 



By using Cauchy Schwartz inequality and since the Brownian motions B and W are independent 
we have: 

E ^Z^' ~ Z^'-^f]h' 

< h'[EUt,,+,,XZ,,,Y,':^^,,E,,^^dZt'^^^^^^^ 

< a'h'Em^^Azt''-'-zZ"-'r] 

which implies for i > 2 

E[\\Z^/ - < a^E[\\Z^^''~' - Z^/-^f]. 

Since < a < 1, this shows that (Z(^'*)i>o is a Cauchy sequence and Zt/ e 2.2(^7^ jrOj ^/jg^g 
L'^{il,J-'f ) denotes the set of square integrable random variables -measurable. We set Z^ :— 

.lim zf''eL\n,Fl). 

Remark 2.3 The superscript {t,x) indicates the dependence of the solution {X,Y,Z) on the 
initial date (t,x). To alleviate notations, we omit the dependence on (t,x) of and 
(Tj^N,t,x ^ ^N,t,x^ yr,iig,fi ifig context is clear. 

Notations: For a real matrix A, \\A\\ is the Frobenius norma defined by \\A\\ = iJ2i j ^1 jY^'^- 
For a vector x, \x\ stands its Euclidean norm defined by |a;| = |a;ip)^/^. 

In the next computations, the constant C denotes a generic constant that may change from line 
to line. It depends on K, T, a, |6(0)|, ||cr(0)||, |/(t,0,0,0)| and 0, 0, 0)||. 



3. Upper bound for the time discretization error 

The following lemma gives an estimation of the second moment of the solution of the BDSDE 
(2.2): 

Lemma 3.1 Under assumptions (H1)-(H2), we have 

E[\Y,-Yt^^,\'']<Ch{l + \x\''). 

Proof: 

From the equation (2.2), we have Vn = 0, .., — 1 and Vs € [in, tn+i) 



E[\Y,-Y,^^,\^] - E 



f{u,X^,Y^,Zu)du 



ZudWu 
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Then 



E[\Ys - Yt, 



< CE 
+ CE 
+ CE 



du 



g{u,Xu,Yu,Zu)dBi^ 



ZudWu 



and so 



E[\Ys~Yt^^,\^] < C E[\f{u,X^,Yu,Zu)\^]du 



+ C E[\\g{u,Xu,Yu,ZuW]dv 



+ C / E[\\Zu\\^]du. 

J s 

From inequalities (3.1)-(2.3)-(2.4) and Assumption (H2), wc deduce that 
E[\Ys~Yt^^,\''] < C f + E[\Xu\^ + \Yu\^ + \\Zu\\''])du 



(3.1) 



< C(h + E[ sup \Xs\^]h + E[ sup \Y,\^]h + E[ [ \\Z,\\^ds]h 

^ 0<s<T 0<s<T Jo 

< Ch(l + \x\^). 



□ 



The following theorem states an upper bound result regarding the time discretization error. 
Theorem 3.1 Assume that the hypothesis (H1)-(H2) holds, define the error 



Error n{Y,Z):= sup E[\Yt - V/" f] + E[ [ \\Z^ - Ztfdt], 

0<t<T Jo 

where Y^ and Z^ are given hy (2.11). Then 



(3.2) 



N-l 



ErrorN{Y,Z) < Ch + Cj^ E\\Zf - Z^J\^ds 



ji=0 

^-1 /•t„+i 



E\\Zs~ Zt^^.W^ds. 



(3.3) 



Proof: We first define the following quantities: 

:= (X3,y,,Z,),Vse [t„,i„+i),Vn = 0,..,iV-l. 
Vse [t„,t„+i),Vn = 0,..,Af-l: 

" jy^^ := Yt - Y/", 8Zf := Z, - Zf , 

Sft ■.= fit,et)-fit^,e^), 

[ Sgt :-g(t,0O-.9(Wi,e^+i). 



(3.4) 
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From the definition of tlie BDSDEs (2.2) and (2.11) we have Vt S [t„,t„+i) 

Jt Jt Jt 

and 

Jt Jt Jt 

Then we get: 



J/sds + / SgsdBs - i 
It Jt Jt 

Using the Generalized Ito's Lemma [26], we obtain: 



Then taking the expectation we find 



At:=E[\SY,''f]+J^ E[\dZ^\']ds~E[\SY,^^J'] ^ 2 E[{SYf ,S^)]ds 



+1 



J^''^*' E[\Sg,\^]ds (3.5) 



From the assumption (H2) 



^ E[\Sgs\^]ds <Kh^+K E[\Xs - X^^, ^ds 

K r'\\Y,-YZj'']ds + a'E[ f^^' \\Z, - E^jZ^JW^ds] . (3.6) 



We use the Young's inequality, for a positive constant e (to be specified later) and then we use the 
projection property of the conditional expectation and Jensen's inequality to obtain 

+ {l + e)E[fj"^ \\E,^^AZs]-E,^^AZZWds\ 

< {1 + -)E[ f"^' \\Zs-Zt^^,\\'ds] 
^ Jt 

+ il + e)E[J^'""'\\Z,~Zi'j\'ds]. 
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Inserting , using the Young's inequality for a positive constant e' (to be specified later) and 
plugging the last inequality in (3.6), we get 

E[\5g,\'']ds < Kh^ + kJ*"^' E[\X, - X^^J'] + K j'"^' E[\Y,, - Y,''^J']ds 
+ (f + i)a2 J'"^" ^[||^^ _ Zt„^,\\^]ds + (1 + e)(l + e')a' J^"*' E[\\SZ^\\^]ds 



1 



+ (l + ,)(l + ^)aM*" E[\\Zf ~ Z,Z\\']ds. (3.7) 



e' 



We set a' := (f + e)(l + e')a^. We choose e and e' such that a' G (0, 1). This is possible since 
G (0, 1). Then we use the inequality 2ab < a? + and the assumption (H2) in (3.5), 

At < Y^J^ E[\SYff]ds + ^^J^ E[\SM']ds + J^ E[\Sgs\']ds 
^ Em:'\']ds + ^^C{h' + E[\X, xZ?]ds 

+ J'"^' E[\Y, - Y,^^^ f]ds + J*"^' E[\\Z, - ZZ\\']ds} + j'j^' E[\5gA^]ds 



+ 1 



EWs-YZj^\ds^2 E[\\SZf\\^]ds + 2 i?[||Zf - 

E[\Sgs\^]ds, (3.8) 



where we inserted Z^ in the third inequality (3.8). 
After that we plug (3.7) in the inequality (3.8): 

4C . 



A < J^J^""E[\6Yj'\']ds+^-^C{h'+J^""E[\X,-XZnds 

E[\Y, - Yf^J']ds + 2 - zZl^ds} + (i-^) E[\\SZ^r]d.s 



+ +Kh' + K E[\X, ~ ^wJ']^^5 + K ^[l^^. - Yt^+f]ds 

+ (1 + -Ja^ J'"^' E[\\Z, - ^t„+JP]ds + a' J^'"*' E[\\SZf\\^]ds 



+ {l + e){l + ^)a' f" E[\\Z^ - Z,Zr]ds. (3.9) 
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We define A'^ hy At- ^ E[\5Zf\'^]d3 



2 Jt 



A[ E[\6Yt^\']+^-^ J^"" E[\6Z^\']ds-Em':j'] 

< C J'"^' E[\SYff]ds + Ch" + C J'"^' E[\X, - X^J^]ds + C j'"^' E[\Y, - Y^^^^jds 

+ C E[\Xs XZ^, \']ds + C J'""^ E[\Y, - Y^^^ \^]ds 

+ C^*"^'i?[||Zf -Z,^|nds + C^*"^'i?[||Z,-Zt„^J|2]ds, (3.10) 

where C is a generic positive constant depending on a' . 
Then, it is known that 

E[\X, - XZ?] < Ch, E[\Xs - XZJ'] < Ch, 
and using the lemma 3.1, we get 

E[\Ys - n < C{E[\Y, - Yt.^. P] + E[\Yt„^, - Y^^l^^ p]} < C{h + E[\Yt^^, - Y^^^ p]} 
This gives 

A[ < C 1'"^' E[\SY^^\^]ds + Ch' + ChE[\SYt'^J'] 



tn + 1 ptn + l 

N ^7^1121^^ I /-i I J7\\\^ ry ||2l 



E^Z^: ^Zl\Y\ds^C E\\\Z,-Z,^^,\\'\ds (3.11) 



From (3.10)-(3.11), we get 

E\\bYl'\^\ < A't + E[\SYt^^J'] 

= E[\6Yt^f]+^-^p'E[\SZ:^f]ds 

< C^"^'s[|,5r/|2]ds + B„, VtG [i„,t„+i), (3.12) 

where we set 

B„ := E[\SYt':^j']+ChE[\SYt';^J']+Ch' 

+ C J^"^' E[\\Zf - zZ\\^]ds + C J^"^' E[\\Zs- Zt,^^,\\^]ds. (3.13) 

Using Gronwall Lemma, we have 

EilSYt'^f] < B„e^^ Vt e [tn,tn+i). (3.14) 
Using now (3.14) in (3.12), we get for h small enough: 

E[m'^\'] + ^ 1^'"" E[\SZf\']ds < il + Che^'^)B^ 

< {i + Ch)Br,, yte[tn,t„+i). (3.15) 
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By taking t = t.^ in the last inequality, we obtain 



Then 



EmZ?\ + ^^ E[\5Z^\^]ds<{l + CK){E[\5Y,^^J^]+ChE[\5YZj^] 

Ch^ + C f E[\\Z^ ~ zZ\\^]ds + C f E[\\Zs~ Zt„^,\\^]ds}. (3.16) 



E[\SY,'^\'] + E[\SZ^\']ds < (1 + Ch){E[\SY,^^J^] 



+1 



+ Ch^ + C I E[\\Z^ ~Z^J\^]ds + C E[\\Z,-Zt^^,\\^]ds}. (3.17) 



Iterating the last inequality and since Nh = T, we obtain Vn G {0, .., — 1} 

1-a' 



^^Il^y*!!!'] + / £^[l'5^f < (1 + Chr{E[\SY^\'] + Ch 

Cj] / -Z/^|p]ds + C^ / £;[||Z,-Z,„,J|». (3.18) 



n=0 n=0 *" 

From equalities (2.2) and (2.8) and using the Lipschitz property of $ (Assumption (H2)-iii), 
wehaveVne {0, ..,A^-1} 



1 — a 



+ 1 



EilSY.'^n + E[\SZf\']ds < (1 + Ch)''{Ch 



/ i?[||Zf -Z,^|p]d. + C^ / E[\\Z.,-Z,,^^,\\']ds}. (3.19) 

„^ — n in. ^ — n ''in. 



Now we sum up the inequality (3.15) for t = tn 



J2E[\6Y,y] + l^ j E[\6Z^\'']ds<{l + Ch){Y,E[\5YZj''] 
ChY,E[\5YZj^ + Ch + cY, / i?[||Zf 

n=0 n=0 

Cj] / E[\\Zs-Zt^^,W']ds]. (3.20) 



Then 



n=0 



1-a' '■^ 



/ S[|(5Zf |2]ds < [(1 + Ch) + (1 + Ch)Ch]E[\5Yr^\^] 
Jo 



2 



N-l 



n=l 

E[\\Z^ -ZZ\\']ds + cY, E[\\Z.,-Z,„^,f]ds, (3.21) 
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which gives for h small enough and using the Lipschitz property of $ (Assumption (H2)-(iii)) 

1 / /-T N-1 N-1 

i-^/ E[\SZf\']ds<C{h + hY^E[\SY,^J'] + Y: / E[\\Zf - zZ\\']ds 

^-1 /•t,.+l 

+ \ E{\\Zs^Zt^^^,\\^\ds). (3.22) 

Moreover summing up (3.19) and using Nh = T 

n=l ri=0 •^*'> 

+ CTJ2 E[\\Zs~Zt„^,\\'']ds}. (3.23) 

n=0 

From inequality (3.22), we deduce that 



f E[\SZ^f]ds < C{h+J2 f""E[\\Z^-ZZ\f]ds 
Jo Jt„ 

J2 / E[\\Zs-Zt„^,f]ds}, (3.24) 



n=0 

On the other hand, using (3.15) we have 

rt,^+l 



E[\dY,^\^] < il + Ch){E[\SY,';^J^] + Ch' + C JJ^^ E[\\Zf ~ Z^J\^]ds 

+ C^"^'s[||Z, -~Zt„^Jp]ds}, Vie [t„,Wi)- (3.25) 



Plugging (3.19) into (3.25) 



E[\SYfn < Ch + CY^I E[\\Z^-Z^X\ds 

n=0 ■ 



Vte [t„,t„+i). (3.26) 



Then by taking the supremum over t € [t„,t„+i) in the last inequality, we obtain 
sup E[\5Y,''\^] < Ch + Cj2 E[\\Zf ^ Z^J\']ds 

t„<t<i„ + i Jt„ 

^-1 /.t„ + i 



E[\\Z,-Zt„^,f]ds. (3.27) 

n=0 -^i" 



Then 



sup < Ch + CY E[\\Z^ - Z^\\'']ds 

0<t<T ^Jt„ 



/ i?[||Z,-Z,„^J|2]d5. (3.28) 



n=0 
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Inequalitites (3.28) and (3.24) give together 

sup E\SYt^f+ [ E\6Z^fds < Ch + Cy [ E\\Z^ - Z^^lfds 

0<t<T Jo ~L Jt„ 



n=0 

+ E\\Zs~ Zt,^_^,fds. (3.29) 



n=0 



□ 



4. Path regularity of the process Z 

The purpose of this section is to prove i^-regularity of the Z component of the solution of the 
BDSDE (1.2) and then to give the rate of convergence of our numerical scheme. For this end, we 
need to introduce the Malliavin derivatives of the solution. This will allow us to provide a repre- 
sentation and a regularity results for Y and Z that will immediately imply the rate of convergence. 
We recall quickly the tools on the Malliavin calculus in the context of BDSDEs introduced in Par- 
doux and Peng [25] . Pardoux and Peng in [26] have skipped details of this part considering that it 
is just a natural extension of the work on standard BSDEs [25]. For the sake of completeness, we 
give some details which are crutial to obtain regularity result of the process Z and we shall give 
all the technical proofs in the Appendix. 

4-1. Generalities 

For any random variable F of the form F = f{W{hi), ...,W{hn), B{ki), B{kp)) with / e 
Cb°°(M"+P,M), hi,...,h„ e L2([o,T],R''),fci,...,fcp e L^{[0,T],W), where 

W{h,) := [ h,{s)dWs, B{k,) / kj{s)tB,. 

Jo Jo 

We let 

D,F := J2 ^..f(w{hi), Wihn); B{ki), B{kp)\ h,{s), 0<s<T, 

where V^/ is the derivative of / with respect to its i-th argument. 

{DsF)s is the Malliavin derivative of F with respect to the Brownian motion W. 

We denote by S, the set of random variables of the above form. For such F, we define its norm as: 

\\F\\,,2 := {E[F^] + E[ / \DsF\^ds]}\ 
Jo 

We define the Sobolev space: 

We introduce the following notations: 

iS^([si,S2],ro^'^) is the set of processes Y = (y„,si < it < S2) such that Y € §^([51,53]), Y^ € ©i-^, 
1 < i < k, si < u < S2 and 

||>"||i,2 -.^{Eif ' \Yu\''du]+E[f ' f ' \\DgYu\\^dude]}i <oo. 

J J S^ J S^ 
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-^fcxd([*i' •^s],®"^'^) is the set of processes Z = (Z^^si < u < S2) such that Z G IHI^^^([si, S2]), 
Zi;^ e ID)i'2 X < i < fc, 1 < j < d, si < w < S2 and 



\\Z\\i,2 ■■={£[[ ' \\Zufdu\+E[r I ' \\DeZ.^fdude]Y^ < cx). 



S2([si,s2],Di.2) :=5|([si,S2],Dl'2)xA^L,([si,S2],Dl'2). 

In order to siniphfy notations, we will omit the dependence on [si, 82] if the time interval is equal 

to [o,r]. 

We define L^([t, T], D^'^) as the set of progressively measurable processes (ws)t<s<r such that : 

(i) v{s,.) e ID)^'2, for a.e. s e [t,T], 

(ii) (s,u)) — > Dv{s,w) e L^i[t,T] x n), 

(iii) E[J^ \vs\^ds] + E[J^ \DuVs\^duds] < 00. 

We denote by L^{[t,T],I])^-^ x D^'^) := L^{[t,T],]n)^'^) x L^{[t,T],]n)^'^). 
4-2. Malliavin calculus on the Forward SDE's 

In this section, we recall some properties on the differentiability in the Malliavin sense. We shall 
assume 

(H4(i)) b e Cl{R'^,R'^) and a G {R'^ ^W^""^) 

(H4(ii)) b e Cb2(M'*,R'*) and a e CU^'^ ,'^'^'"^) 

Then under (H4(i)), Nualart [23] stated that X*'^ e D^-^ for any s e [t,T] and for ^ < A: the 
derivative is given by: 

(i) DlXl''' = 0, for s<r <T, 

(ii) For any t < r < T, a. version of r < s < T} is the unique solution of the linear SDE 

Dlxl-^ = + r Vb{Xl;^)D[xl;^du + E Va'(X*'-)7^[X*'-dI^^, 

where {cr'^)i=i^...^k (resp. {cr^)i=i^...^k) denotes the i-th (resp.l-th) column of the matrix a. 
Moreover, if Assumption (H4(ii)) is satisfied, then, from Nualart [23], D\Xl'^ e D^'^ for all 
r,s < T. For all w < T and I' < k, we have 

dIdIxY-' = if s < i)Vr, 

and for all s > V r a version of D\, D\.Xl'^ is the unique solution of the SDE: 

d 

4=1 



[Y,y{{yby{xr))Dixr{Dixry + ^b{xl;^)DiDixl; 

E/' [Y.^{^^\Xli-)yDixl:-{D[Xl;-y +ya\Xl;-)DiD[Xl;^ 



i=i - ' j=i 
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where ((V6)-^)j=i_..._fc (resp.((V(T'(X^'^))-')j^i ^j) denotes the j-th column of the matrix (V6) 
(resp. (Va'(X*'^))) and ((Li^X*'^)^)j=i,...,fe denotes the j-th component of the vector (Li^X*>=^). 
The following inequalities will be useful later. For the proofs, we refer to Nualart [23] for example. 
From Lemma 2.7 in [23] applied to X and D^X and any 0<r<s<T, we have the following 
inequalities 

E 



sup ll-Ds^uj 

0<ti<T 



E 



sup I 

s\/r<u<T 



<C(l + |x|^), (4.1) 
<C|s~r|(l + |x|P), (4.2) 



and the similar argument for DrDgX shows that there exists C > such that 

e\ sup WDrDsXuW"] <C{l + \x\^P). (4.3) 

4-3. Malliavin calculus for the solution of BDSDE's 

Now, our aim is to study the differentiability in the Malliavin sense of the solution of the BDSDE 
(2.2). We start with the following lemma which shows that a backward Ito integral is differentiable 
in the Malliavin sense if and only if its integrand is so. We recall that Pardoux and Peng [25] 
proved that the result holds for the classical Ito integral. 

Lemma 4.1 Let U G H^([t,T]) and I,{U) = UrdW^J = Then, for each 9 e [0,r] we 

have Ug G D^'^. if and only if Ii{U) £ D^-^ji = 1, ...,d and for all 9 G [0,r], we have 



DgI,{U) = / DgUrdW^ + Ug, 9>t, 

Je 

Dgh{U) = DgUrdW^, 9 < t. 



For backward Ito integral, and since the Malliavin derivative is with respect to the Brownian 
motion W, we have the following result : 

Lemma 4.2 Let U £ ml{[t,T]) and Ii{U) = UrdBl.,i = 1,...,L Then for each 9 e [0,T] we 
have Ug € D^'^ if and only if Li{U) S D^'^, i = 1, I and for all 9 G [0, T], we have 



DgL,{U) = / DgUrdBi, e>t, 

DgL,{U) = ^ DgUrdBl, 9<t. 



For later use, we need to prove the a priori estimates for the solution of the BDSDE which is the 
object of the following proposition. 

Proposition 4.1 Let {(1)^ , f^ , g^) and (0^, .g^) be two standard parameters of the BDSDE (2.2) 
and (Y^,Z^) and (y^,Z^) the associated solutions. Assume that Assumption (H2) holds. For s G 
%T],putSY, ■.^Y}-Yl52fs := f\s,X,,Y^,Zl)-f\s,X,,Y^,Zl) andS^gs := g\s, X,,Yl Z^y 
g'^{s,Xs,Y^,Zl). Then, we have 

\m\liiit,T]) + \\m\li^^ilt,T])<CE[\5YT^ (4.4) 
where C is a positive constant depending only on K , T and a. 
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Proof: The proof of this result is classical in the BSDE's theory so we omit it. □ 

Now, we study the differentiability in the Malliavin sense of the solution of the BDSDE which 
is technical. To our knowledge, it does not exist in the literature. We have to precise that Pardoux 
and Peng [26] have skipped details considering that it was just an easy extension of the work 
on standard BSDEs [25]. We show that the derivative is a solution of a linear BDSDE as Peng 
and Pardoux [25]have given for the standard BSDE's, see also El Karoui Peng and Quenez ([9], 
Proposition 5.3)). 

Proposition 4.2 Assume that (H2)-(H3(i))-(H4(i)) hold. For any t e [0,T] and x e W' , let 

{{Ys,Zg),t < s < T} denotes the unique solution of the BDSDE: 

= 0(X^")+ / f{r,Xl^^,Yr,Zr)dr+ f g{r, X^^^ ,Yr, Zr)tBr - f ZrdWr, t<s<T. 

Us Us Us 

Then, {Y,Z) e B^{[t,T],D^-^) and {DgY,, DgZ,;t < s,6 < T} is given by: 
(i) DgYs = 0, DgZs = for all t < s < 6 <T 

(a) for any fixed € [t,T], < s < T and 1 < i < d, a version of (DgYg, DgZg) is the unique 
solution of the BDSDE: 

DlYs - V$(X^")i?^X^" + ^ {yj{r,Xl^-,Yr,Zr)DlXl^-yr 

/T d 
hyfir, Xl'\ Yr, Zr)DlYr + J2 ^^^fir. K'^ Yr, Zr)DlZl)dr 

+ {^.g''{r,Xl'\Yr,Zr)DlXl.^- + Vyg'\r,Xl^\Yr,Zr)DlYr)dB^ 

n=l •''^ 

I «T d ^ „T d 

+ E/ Y.i^-'9''{r,Xl-\Yr,Zr)Dlzi)dB': - Y^^oZldW}. (4.5) 

where {z^)i<j<d denotes the j-th column of the matrix z, {g"')i<n<i denotes the n-th column of the 
matrix g and B = {B^ ^ i?'). 

Proof: See Appendix. □ 

The second order differentiability in the Malliavin sense of the solution of the BDSDE will be 
given in Appendix. 

4.4- Representation results for BDSDEs 

In this subsection, we will prove a representation result of {Z, DZ) which will be useful to prove 
the rate of convergence of our numerical scheme. 

Proposition 4.3 Assume that (H1)-(H4) hold. Then : For t < s < T, we have 

DsYs = Z„ (4.6) 
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and 

<^^(i+Nn- (4-7) 

For li,l2 ^ d, t < s < T , we have 

Di-D['Ys=D\-Zi\ (4.8) 

and 

ll^i^^ll|,,([t,T]) <C^(1 + NI')- (4-9) 

Proof. To simplify the notations, we restrict ourselves to the case k = d = 1. 
1. Notice that for t < s 



dW, 



where := (X^^^, y^, Z^). 

It follows from Lemma 4.1 and Lemma 4.2 that, ior t < 9 < s 

DgYs = Zg- f (w,.J{r,i:r)DeXr + Wyf{r,j:r)DeYr + WJ{r,i:r)DgZr)dr 

'J 6 

- f hx9ir,^r)DeXr + Vyg{r,j:r)DgYr + V,g{r,j:r)DgZr)^r + I DgZrdWr 
JO J 6 

Then by taking 6 = s, it follows that equality (4.6) holds. From (7.1), we deduce that (4.7) holds. 
2. Notice that for 6* < t < s 

DgY, - ^"^'"Z {y-f{r,^r)DgXr + ^yf{r,i:r)DeYr + ^,f{r,^r)DgZr)dr 

- Z {yx9{r,^r)DeXr + Vyg{r,^r)DgYr + V,g{r,i:r)DgZr)'iiBr + DgZrdWr. 
It follows from Lemma 4.1 and Lemma 4.2 that, for 9 < t < v < s 

D^DgYs = DgZ,- [ D4i:r)*[Hf]{r,Er)Dg{^r)dr- f Vf{r,^r)DM^r)dr 

J V J V 

- j D.,{^rT[Hg]{r,^r)Dg{^r)tBr ~ j Vg(r, S,)Z?,i^e (S,)^ 

J V J V 



+ / D.DgZrdWr. 
J V 

Then by taking w = s and t = 9, it follows that equality (4.8) holds. We have from estimate (2.4) 
and inequality (4.3), that for each v <T and 9 < T 

E[ sup \D,DgYs\''] +E[[ \D,DgZs\^ds] < C(l + \x\*). (4.10) 

t<s<T Jt 

and then by taking v — s and t — 9 wc deduce that (4.9) holds. □ 
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4-5. Path regularity 



In this subsection, wc extend the resuh of Zhang [29] which concerns the L^-regularity of the 
martingale integrand Z . Such result is crucial to derive the rate of convergence of our numerical 
scheme. We start with the following proposition which gives an upper bound for 



E 



sup \Yr-Ys\ 



and E 



, t<s<u<T. 



Proposition 4.4 Assume that (H1)-(H4) hold. Then for t < s < u < T , we have 

E sup < C{l + \x\'^)\u~ s\, 

r^[s,u] 



E 



< Cil + \x\'')\u-s\ 



(4.11) 
(4.12) 



Proof. To simplify the notations, we restrict ourselves to the case k = d = I = 1. 

(i) Plugging inequality (4.7) in the estimate (2.5), the result (4.11) holds. 

(ii) From Proposition 4.3, we have 

\Zu-Z,\^] < CE[\DuYu-DsY,\^] + CE[\DsY,-DsYs\% 
From the definition of the BDSDE (4.5), we have 

DuYu ~ DsYu = VHXt){DuXt - DsXt) + I (Vxfir, I],)(D„X, - DsXr) 

J U 

+ I hyf[r,^r){DuYr-D,Yr)+y,f[r,^r){DuZr-DtZr))dr 



(4.13) 



dr 



V,g{r, llr){DuXr - D^Xr) + Vyg{r, S,)(i5„y, - DsYr)jdBr 

f h,g(r, ^r){DuZr - D,Zr))tRr - f {DuZr - D,Zr)dWr. 



A. Bachouch, M.A. Ben Lasmar, A. Matoussi, M. Mnif/ 20 

Applying the generalized Ito's formula, we obtain 

- 2 / V,/(r, ^r){DuXr - D,Xr){DuYr - D,Yr)dr ~2 f Vyfir, ^r){DuYr - DsYrfdr 

J u J u 

- 2/" V,f{r,i:r){DuZr-D^Zr){D^Yr-D,Yr)dr 

J U 

- 2 j \7,g{r, - DsXr){DuYr - D,Yr)tu3r 

J U 

-2/" \7yg{r,j:r){DuYr-DsYrftlBr 

J u 

- 2 / V,5(r, ^r){DuZr - DsZr){DuYr - DsYr)tlBr 

J u 

+ 2 ( [DuZr-D,Zr){D^Yr-DsYr)dWr 
J u 

J U 

+ / \D^Zr- DsZr\^dr 

J u 

From inequalities (7.1) and (4.1), using the Burkholder-Davis-Gundy's inequality and Assumption 
(H2), the stochastic integrals which appear in the last equation disappear when we take the 
expectation. By Young inequality, we obtain, for e' > 

E[\DuY^ - DsYu\^] +E[( \D^Zr - D,Zr\^]dr < E[\V^Xt){DuXt - DsXt)\^] 

J U 

+ 2E[f V^f{r,i:r){DuXr~ DsXr){DuYr- DsYr)dr\ 
J u 

+ 2E[( \/yf{r,^r){DuYr~ DsYrfdr] 

J u 

+ 2E[f VJ{r,i:r){DuZr-D.uZr){DuYr^D,Yr)dr] 

J u 

+ C{l + \)E[j V,5(r,E,)'|i?nX.- As^rl'dr] 
e Ju 

+ C{l + ^)E[( ^yg{r,^rf\DuYr- DsYr\^dr] 

^ J U 

+ {l + e')E[( \/,g{r,j:rf\DuZr-DsZr\^dr]. 

J U 



A. Bachouch, M.A. Ben Lasmar, A. Matoussi, M. Mnif/ 



21 



Hence by using Assumption (H2) and Young inequality, we have for e, e' > and C > 0, 

E[\DuYu - DsYu^"] +E[f \DuZr - DsZr\^dr] < K^E[\DuXt - DsXt^"] 

J U 

+ 2KE[ / \DuXr - D,Xr\'^dr\ + AKE[ / \DuYr - D^Y^^^dr] 

J U J U 

+ KeE[[ \D„Yr-DsYr\^dr] + —E[[ \D^Zr - D^Zr^dr] 

J U ^ J U 

+ CK\l + ^)E[f \DuXr-DsXr\^dr]+CK\l + ^)E[f iD^Yr - D^Yrl^dr] 

^ J U ^ J U 

+ {l + e')a^E[f \DuZr-DsZr\^dr]. 

J u 

Then, we obtain 

E[\DuYu - D.Y^f] +E[[ \DuZr - DsZrfdr] < K^E[\DuXt - DsXrf] 

J U 

+ K{2 + KC{l + -^))E[( \DuXr-DsXr\^dr] 

+ {K^C{l + -, ) + {A + e)K)E[l - i^.F.pdr] 

^ J U 

+ ((1 + e')a^ + —)E[ I \DuZr - D^Zrl^dr]. 
e Ju 

For e large enough and e' small enough, we have (1 + e')a^ + < 1. From inequality (4.2), we 
deduce that 

E[\DuY,, - DsYul^] < c({l + \x\^)\u^s\+E[f \DuYr ~ D^Yrl^dr]) , 

J U 

where C is a positive constant. From Gronwall's lemma we have 

E[\D^Y^ - DsY^\^] < C{1 + \x\-')\u - s\. (4.14) 

Since {DsYu)s<u<t satisfies the BDSDE (4.5), inequalities (2.5)-(4.7) hold for {DgYu, DsZu)s<u<t 
and yield 

E[\D,Yu - DsYs^"] < C(l + \x\^)\u - s\. (4.15) 
Plugging (4.14) and (4.15) into (4.13), we obtain (4.12). □ 

4-6. Rate of convergence for the BDSDE 

The following theorem shows the rate of convergence of our numerical scheme (2.9)-(2.10). 

Theorem 4.1 Under Assumptions (HI) and (H2), there exists a positive constant C (depending 
only on T, K, a, |6(0)|, ||cr(0)||, |/(t,0,0,0)| and \\g{t,Q,Q,Q)\\) such that 

ErrorN{Y,Z)<^{l + \x\''). (4.16) 
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Proof: We consider a density function qi (resp. (72) defined on (resp. 
compact support. For to S N, we define 

(&",CT")(a;) = m'^ j {b,cr){x)qi{m{x - x))dx, 



X R 



kxd\ 



with 



d+k+kxd 



{f , 9){v)q2{m{v - v))dv, 



where v = {x, y, z) and v = (x, y, z). For m large enough, the functions 6™, ct™, /™ and are K- 
Lipschitz and C^. Also (6™, ct™, /™, 5™) converge uniformly to (6, tr, /, g) on every compact and so 
for TO large enough, we have |6™(0)-fo(0)| < 1, | |ct™(0) - cr(0)| | < 1, 0, 0, 0) - /(i, 0, 0, 0)| < 

1 and ||5™(t,0,0,0) - 5(^,0, 0,0)|| < 1. Besides, we note that (6™, cr™, g™) converge in 
to (6,CT,/,.g). We denote by the solution of (2.1)-(2.2) with {b,a,f,g) is replaced 

by (6™, (T™, /'", g™). From Lemma 7.1 and inequality (4.4), we have (X'",F™,Z™) converge to 
(X, y, Z) in §^ X Sfc X From Proposition 4.4, we have 



E 



\ZT - Z^:f] < C(l + \x\^)\s -U\<^il + Ixp), 



for all s e [U, U+i] and < i < iV - 1 where C depends only on T, K, 6"(0), cr'"(0), /""(t, 0, 0, 0) 
and .g™(t, 0, 0, 0). Integrating over s and summing over i, it yields 

^-1 - rU+i 



i=0 

Sending m to infinity, we deduce that 

^-1 nU + i 



i=0 



\Zs — Zt \\ ds 



<§(i + N^)- 



Now applying the last inequality for Z , we obtain that 

N-1 



N-1 

ml 



\Z^-Z^\\^ds 



Similarly, we have 



i=0 



N-1 



<^(i + N^). 



■1=0 



'ds 



Plugging (4.19) and (4.18) in (3.29), we deduce inequality (4.16). 



(4.17) 



(4.18) 



(4.19) 



□ 



5. Numercical scheme for the weak solution of the SPDE 

Most numerical work on SPDEs has concentrated on the Euler finite-difference scheme (see [13], 
[14] , [12]), on finite element method (see [28]) and also on spectral Galerkin methods (see [15] 
and the references therein). Here, we follow a probabilistic method based on the Feynman-Kac's 
formula for the weak solution of the semilinear SPDE's (1.1) based on BSDE's approach (see [5], 
[22]). We consider a weak Sobolev solution of such SPDE in the sense that u shall be considered as a 
predictable process in some first order Sobolev space. Therefore, we shall improve the convergence 
and the rate of convergence of the L^-norm error of such solution by using the convergence results 
on BDSDEs proved in section 4 and an equivalence norm result given in Barles and Lcsigne [4] 
and Bally and Matoussi [5]. 



A. Bachouch, M.A. Ben Lasmar, A. Matoussi, M. Mnif/ 



23 



5.1. Weak solution for SPDE 

Since we work on the hole space W^, we need to introduce a weight function which is integrable 
and satisfies /^^(l + \x\'^) p{x)dx < oo. For example, we can take p{x) = or p(x) — e^'^'). We 
assume the following integrability conditions on ($,/,(;): 

(Hp) (i) / \<^(x)fp{x)dx<^, 

(ii) / / |/(t,a;,0,0)|V(a;)da;di < oo, 

(iii) \ [ \g{t, x,0, 0)1'^ p{x)dxdt < oo. 
Jo Jm^ 

We denote by L'^(M.'^,dx) the basic Hilbert space and we employ the usual notation for its scalar 
product and its norm {u,v) — J^^ u{x)v{x)dx and ||w||2 = {u,u)^. Then, we define by H;^(M'*) 
the first order Dirichlet space and its norm ||m||hi(R'') = (II""!!! + ll^^'^lll)^- We define also V :— 
C^{[0,T])®C'^(R'^) the space of test functions where C;f([0, T]) denotes the space of all real valued 
infinite differentiable functions with compact support in [0,r] and C^(R'^) the set of C^-functions 
with compact support in M''. 

We introduce Ht the space of predictable processes {ut)t>o with valued in H;^(M'') such that 











sup ||wt|l2 


+ E 


[ \\Vuta\\ldt 


1 < OO 


'-o<t<r J 









Definition 5.1 We say that u G "Ht is a weak solution of the equation (1.1) associated with the 
terminal condition $ and the coefficients {f,g), if the following relation holds almost surely, for 
each if £ V 



.), dsifis, .))ds + 5(u(s, .), ^{s, .))ds + {u{t, .), ^(t, .)) - (*(■), ^{T, .)) (5.1) 
(/(s, ., u{s, .), {Vua){s, .)), <f{s, .))ds + Y, {gis, ; u{s, .), (Vwa)(s, .)), (^(s, .))dBl, 



where 



£{u,ip) ^ {Lu,ip) ^ I {{\/ua){\/fa) + Lp\/{{^a*\/(T + b)u)){x)dx 



is the energy of the system associated with the SPDE. 
From Bally and Matoussi [5], we have the following result: 

Theorem 5.1 Under (H1),(H2), (H4) and (Hp), there exists a unique weak solution u € Ht 
of the SPDE (1.1) associated with the terminal condition $. Moreover, u{t,x) — Y^'^ and Zl'^ — 
Vuta, dt(i)dx®dP a.e. where Zl''^)t<s<T is the solution of the BDSDE (1.2). 



5.2. Numerical Scheme for SPDE 

Let us first recall that (X^, y^, Z^) denotes the numerical Euler scheme of the FBDSDE's (1.2) 
given in (2.7)-(2.8)-(2.9)-(2.10). The numerical approximation of the SPDE (1.1) will be presented 
in the following lemma: 
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Lemma 5.1 Let x (^M!^ and t,tn S tt such that t < tn. Define 

ul {x) Yll''-^ and < (x) (5.2) 
Then (resp. v^) is J^^ rp-measurable and we have 

Proof: From the Markov property of and , the random variables and tj^ arc J-f j, 
measurable. From the definition of and u^, wc have 

"t„ (^t„ ) = and (X^; ) = 

From (2.9), (2.10) and by taking {t,x) = {tn,Xl'^), we obtain : 

From the flow property, we have X^^'*"'"^*" = X/^'*'^, then we obtain 

= E,^ k::;*;'^- AM/; + ff(t„+i, xi!;?, r.^:;*:'^- , [zj ])as„ am/„i. 

Then from the uniqueness of the solution of (2.9)-(2.10) we obtain the result. □ 
5.3. Rate of convergence for the weak solution of SPDEs 

We give a norm equivalence result which was already proved by Barles and Lesigne [4] and Bally 
and Matoussi [5] when b € (^^^(IR'', M'') and a e Cl{W^ ,W^^'^). Here, we do not assume that a 
satisfies the uniform ellipticity condition. 

Proposition 5.1 Under (H4)-(Hp), there exist two positive constants Ci and C2 such that for 
every t < s <T and (j) G L^(M'' x Vis, p{x)dx ® dPs), we have 

Ci f EBmx)\]p{x)dx< f EMXm\]p{x)dx<C2 I EBmx)Mx)dx, (5.3) 

JR^i JK'' JWL<i 

where Eb is the expectation under Pb- Moreover, for every ^ G L^{W^ x (0,7") x Qb, p{x)dx <E) 
dt ® dPB) 

Ci [ I EB[\^{s,x)\]dsp{x)dx < I I E[\^{s,Xl''')\]dsp{x)dx 

< C2 [ [ EB[\'i'{x)\]dsp{x)dx. (5.4) 
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We recall that x) = y/'^ and v{t, x) = Z*'^ dt® dx® dP a.e. Wc define the process (u^, vf) 
as follows: 

(x) := '^'^ and (z) := Zf■«■^Vs € [t„,t„+i). (5.5) 
Using (2.11) and following the proof of Lemma 5.1, we obtain 

(X*'-) = -^^^ and (X*'-) = Zf'*-^Vi < e [t„,t„+i). (5.6) 

We define the error between the solution of the SPDE and the numerical scheme as follows: 
Error Miu,v) := sup Eb[ \u^ (x) ~ u{s, x)\'^ p{x)dx] 

0<S<T jRd 

+ H^b[ / \\v^{x)-v{s,x)fdsp{x)dx\. (5.7) 

The following theorem shows the convergence of the numerical scheme (5.2) of the solution of the 
SPDE (1.1). 

Theorem 5.2 Assume that (H1),(H2), (H4) and (Hp) hold. Then, the error Error n{u,v) con- 
verges to as N oo and there exists a positive constant C (depending only on T, K, a, |6(0)|, 
||cr(0)||, |/(t,0,0,0)| and \\g{t, 0,0,0)\\) such that 

(J 

Error m{u,v) < —. (5-8) 

Proof: We take t — to. From the norm equivalence result (see inequality (5.3)), 
Vs G [tn,tn+i) such that s > t, we have 

Eb[ f \u^{x) - u{s,x)\'p(x)dx] < CE[ f \u^{Xln " u{s, X^^^l^ p{x)dx], 

where C is positive generic constant. From equation (5.6), we get 

EB[f \u^{x)-u{s,x)\'p{x)dx]<C f E[\Yf''^^~Y:^^\M^)dx. 

Therefore Theorem 4.1 implies that 

sup Eeif \u^{x)-u{s,x)\^p{x)dx]<^ f (l + \x\')p{x)dx<^. (5.9) 

0<s<T jRd JV Jjjd iV 

From the norm equivalence result (see inequality (5.4)), we have 

N-l ,t 



^Eeif"*' f \\v^{x)-v{s,x)\\^p{x)dxds] 



n=0 



^ ^E^[/ / \w:ixin-v{s,xinrpi^)dxds]. 

n=0 -^K" 

From equation (5.6), we get 



Y,E[f f \\v^{xin-yis,xin\\''p{x)dxds] 

Y,E[ / \\Zf''^^-Zr\\'p{x)dxds], 
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and so from Theorem 4.1 we deduce that 

Y^EbiT'^ I Wv^ix) - vis, x)\\' p{x)dxds]< ^ [ (l + \x\^)pix)dx<^. (5.10) 

From incquahtics (5.9) and (5.10), we deduce that (5.8) holds. □ 
Remark 5.1 Gyongy and Krylov [12] considered the following linear SPDE on [0,T] x M.'^, 

oo 

du{t, x) = (£iu(i, x) + f{t, x))dt + ^^(£2,i'u(i, x) + g{t, x)i)dwl 

i=l 

1 ^ / \ / N d'^u(t,x) du{t,x) ,/,/ \ \oo 

where Ciu[t,x) = ^ a[t,x)ig——- , £,2.iu(t,x) = ^biq[t,x)^- ,1 < i < oo and [b[t,x)^^g)i^Q G 

q.l=l OXlClXq dXq 

e'^, (t,x) e [0,T] X M'^, 1 < g < rf. They approximate the SPDE by 

OO 

du^{t,x) ^ {dlu'\t,x) + f{t,x))dt + Y^{C'i^,u^\t,x) + g{t,x)^)dwl, 

i=l 

Ci. are the approximation of Ci, €2.1 by using finite difference scheme on the space gridGh- 
Their results revolve to prove the existence of the random process u^^^ {t, x), j = l..k for some k > 
s.t. 

u^{t, x) = u^^'^t, x)+J2 ^7«^'^ 2;) + R'^it, x), 

where IS the solution of the SPDE. They assumed that the SPDE is non degenerate and for 
m > fc + 1 + the coefficients are m-times continuously differentiable in x. When they used a 
symmetric finite difference scheme and d = 2, the L"^ -error is proportional to where h is the 
discretization step in space and by the Richardson acceleration, the error is proportional to h"^. 
Compared to their work, our scheme is more general. It converges in the non linear case. Our 
convergence is of order \fh where h is the discretization step in time. However, in our work, the 
rate of convergence does not depend on the space dimension d. 

Remark 5.2 If we assume more regularity conditions on the coefficients and the final condition 
as in Pardoux and Peng [26], namely, $ G C^^ (j^d ^ j^fc ^ ^ j g C^{[0,T] x R'^ x M*-^ x M'^><'^', R'^) and 
gG Cj3([0,T]xR''xR'^xR'*^'^,R'''^'). // (F^*-^, Z*'^)t<,<T is the solution the BDSDE (1.2). Then, 
ut{x) = Yf'^ , y{t,x) G [0,r] X is the unique classical solution of the SPDE (1.1) in the integral 
sense (see [26]). Therefore, we can use (4-16) instead of (5.8) and obtain a stronger result. In fact, 
using (3.2), the estimation on the error (5.1) obtained in the previous theorem is replaced by: 

sup \uf{x)'u{t,x)\'^]+E[ \\v^{x)-v{t,x)\\'^dt\<ch. 

0<t<T Jq 



This last equation gives an estimation which holds for all x G M'' and which is not only almost sure 
anymore. 
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6. Implementation and numerical tests 

In this part, wc are interested in implementing our numerical scheme. Our aim is only to test 
statically the convergence of this scheme. Further analysis of the convergence of the used method 
and of the error bounds will be accomplished in a future work. 



6.1. Notations and algorithm. 



Wc use a path-dependent algorithm; for every fixed path of the brownian motion B, we approxi- 
mate by a regression method the solution of the associated PDE. Then, we replace the conditional 
expectations which appear in (6.1) and (6.2) by L^(r2,P) projections on the function basis ap- 
proximating L^(J7, J-"t^). We compute Z^^ by using I Picard iterations as it is given in an implicit 
way, and then we compute in an explicit manner. Actually, we proceed as in [11], except that 
in our case the solutions and Z^^ arc measurable functions of {X^^^, {ABi)n<i<N-i)- So, each 
solution given by our algorithm depends on the fixed path of B. 



6.1.1. Numerical scheme 



For each fixed path of B, the solution of (2.1)-(2.2) is approximated by {Y^ ,Z^) defined by the 
following algorithm, given in the multi-dimcnsionnal case: 
For < n < iV- 1: 
Vji e {i,..,fc}, 

I 



= Eu 



Vn e {l,..,fc} and V72 e {i,..,d} 



(6.2) 



We stress that at each discretization time, the solution of the algorithm depends on the fixed path 
of the brownian motion B. 



6.1.2. Vector spaces of functions 

At every we select k{d + 1) deterministic functions bases (jpi,n{-))i<i<k{d+i) and we look for 
approximations of Y^ and Z^^ which will be denoted respectively by and , in the vector space 
spanned by the basis {Pji,n{-))i<ji<k (respectively {Vji ,j2.n{-))i<ji<k,i<j2<d) ■ Each basis Pi,„(.) is 
considered as a vector of functions of dimension Li^n. In other words, Pi,n{-) ~ {c(-Pi,n{-),ct € 

As an example, we cite the hypercube basis (HC) used in [11]. In this case, J5i,n(.) does not depend 
nor on i neither on n and its dimension is simply denoted by L. A domain D C M.'^ centered on 
Xf) = x, that is D = Y\i=i{^i ~ a,Xi + a], can be partitionncd on small hypercubes of edge 6. Then, 
D^Ui, . Ai .. where Di^^,,^i^ = (x, - a + ii6, - a + iid] x ... x {x^ - a + 1^(5, - a + idS]. 
Finally we define Pi.n{-) as the indicator functions of this set of hypercubes. 
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6.1.3. Monte Carlo simulations 



To compute the projcetion coeffieieiits a, we will use M independent Monte Carlo simulations of 



andAVKi which will be respectively denoted by X^'"'"' and ^!AV^\m = l, ..,AI 



6.1.4. Description of the algorithm 

Initialization: For n ~ N, take (y^'™) — 
Iteration: For n ~ N — 1, 0: 

• We use / Picard iterations to obtain an approximation of Zt^ : 

• For * = 0, Vji e {1, .., k} and j2 e {1, .., d}, a^^;° = 0. 

• For i = 1, .., /: We compute first -Et„+i i^t!,] appearing in (6.2): 

V?i e {i,..,fe},Vj2 e 



/M,i-1 



1 ^ I 

■ r \ ^ MA—1 m I m 

1 = argmf — ^ - « -Pji 



then we set (•) = ("j^Ln+i-^'ji J2,n+i(-))> ii & {1, -fc}, :'2 e {1, -d}- 

After that, we approximate (6.2) by computing 



m— 1 



.J2 



/ 

EfvN.m N.M/vN,m\ i lM.i-1 m \\ 



Then we set z^f^f^{.) = {af^j^^ri-Pju^A-))^ h e {1, ..fc}, 32 e {1, ..d}. 
• We first compute Et^^^lZ^] appearing in (6.1): 



a 



1 I 

1 = arginf — ^ \afJ^,n-Pn.j2.n " 



m— 1 



then we set 4+ijij2(-) = ("jt!i2,"+i-?'ji J2,"+i(-)), ii € {1, -fc}- ^2 e {1, -d}. 

After that we approximate (6.1) by calculating ctf^ n-- ^^r every ji e {I, .., k}, as the minimizer of: 

M 
m— 1 

/ 

(vN.m N.M,vN,m\ lN,M,I/vN.rn\ a ^ m 



Finally, we define y^'^^(.) as: 



y^f(-) = «n-P..,n(-)),Vjie{I,..,fc}. 



A. Bachouch, M.A. Ben Lasmar, A. Matoussi, M. Mnif/ 



29 



6.2. One dimensionnal case (Case when d — k — I — 1) 
6.2.1. Function bases 

We use the basis (HC) defined above. So we set: 

di = minX™, d2 = niaxX™, and L = -^—z — - 



where S is the edge of the hypercubcs (Dj)i<j<L defined by Dj 
We take at each time t„ 



d + {j-i)s,d + jsyyj. 



^DjiX^j'^) ~ l[d+{j-i)S4+jS){X^j'^),j — 1, ..,L 



and 



Card{Dj) denotes tlie number of simulations of X^^ which are in our cube Dj. 



This system is orthonormal with respect to the empirical scalar product defined by 

M 

M 



m—l 

In this case, the solutions of our least squares problems are given by: 
1 

\ " fvN,m\, M,i-1 m \ 

TJ 2^ Pl.n+l(^t„+j(«jij2,n-Pjij2,n) 



lM,i-l _ 
^l,Ti+l — 



M 

m—l 



M.: 



j^l^piA^t^ )[yn+i (^t„+J— ^ 

m—l 



M 

J- fvN,m\^ N.MivN,m\ , uc(vN,m N,MfvN,m\ N.M.I ivN,m\ 

J^l^PoA^tJ )[yn+i i^tj+r) + h.f[XtJ ,y„^Li (Xt^;J,z„ ' {XtJ ) 

m—l 

9[Xt„^,,yr^u (^t„;j,Wi (^t„;jJAB„ j 



Remark 6.1 FFe note that for each value of M , N and S, we launch the algorithm 50 times and 
we denote by (5^Q''^^'^'*'^)i<i<5o the set of collected values. Then we calculate the empirical mean 
yO,x,N,M empirical standard deviation a^'^^ defined by: 



50 



i=l 



^ 50 

JtEI<;""'"-^o""'"P- (6-3) 



49 . ^ 



Let us finally stress that (6.3) gives us an approximation ofu(0,x) the solution of the SPDE (1.1) 
for t = 0. 
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6.2.2. Case when f and g are linear in y and independent of z 
dXt = Xtifidt + adWt), 

$(a;) = (a- - A'i)+ - 2{x - K2)+ , f{y) = aoy, g{y) = boy 

and we set Ki = 95, K2 = 105, r = 0.01, R = 0.06, Xq = 100, fi = 0.05, a = 0.2, T = 0.25, di = 60 
, (^2 = 200 , oq and bo are fixed constants. 

Let Yexpiicit be the solution of our BDSDE in this particular case. By an integration by parts 
formula we get 

At t=0, we have 

Q X N Ad 

Then, we define Yq ' ' as the numerical approximation of the solution of the BDSDE in this 
case (computed by our algorithm), a^'^^ its standard deviation. In the other hand, we compute the 
solution ^"e^pijcit ^^^^ linear case thanks to the Black and Scholes model, by using the explicit 
formula of the expectation of our Payoff: 



^ e(-o-^fc^)T+6oi3T|^gP^(^^^^.J _ KiN{d2,K,) - 2[xem{di.K,) - K2N{d2.K,)]}, 

where N{.) denotes the standard Gaussian cumulative distribution function and 

I di,« = ^{lnf + (/.+ ^)T}, 
[ d2.R = di^R - aVr, 
for a real non negative constant R. 

For ao = 0.5, b = 0.5 and 6=1 



N=20, Yg's = 2.789 
For ao = 0.5, b = 0.5 and 6 = 0.5 

N=30, Y°'g = 2.869 



M 




1 ' BS ~' I 


100 


2.978(0.365) 


0.067 


1000 


2.792(0.128) 


7 10-4 


10* 


2.807(0.052) 


0.006 




M 




1 *^BS ~^ 1 
^ BS 


100 


2.988(0.446) 


0.041 


1000 


2.855(0.143) 


0.004 


10^ 


2.890(0.053) 


0.007 



A. Bachouch, M.A. Ben Lasmar, A. Matoussi, M. Mnif/ 



31 



6.2.3. Comparison of numerical approximations of the solutions of the FBDSDE and the FBSDE 
: the general case 

Now we take 

$(x) = (x - Kl)+ - 2{x - K2)+, 
f{t, X, y, z) = -6*2 - ry + {y - ^)-{R- r) 
gi{t, X, y, z) — O.lz + 0.5?/ + log{x) 

and wc set B ^ {^i - r)/a, Ki = 95, K2 = 105, Xq = 100, n = 0.05, a = 0.2, r = 0.01, R = 0.06, 
51 = 1,N = 20 , T = 0.25 and we fix di = 60 and ^2 = 200 as in [10]. In this case, assumptions 
(H1)-(H4) are satisfied. 

We compare the numerical solution of our BDSDE (noted again Y ^ ' ' ) and the BSDE's one 
(noted here by Y ^ bsde ): without g and B. 

We finally note before starting the numerical examples that, for the contraction constant taken in 
the following (a = 0.1), our algorithm converges after at most three Picard iterations. 
When t is close to maturity 



M 


^ tio,BSDE\" I 


^ tio y'^ ' ) 


128 


2.925(0.35) 


3.851(0.404) 


512 


2.814(0.178) 


3.711(0.206) 


2048 


2.795(0.094) 


3.684(0.105) 


8192 


2.8(0.054) 


3.690(0.059) 


32768 


2.81(0.025) 


3.699(0.028) 




M 


^ tir^,BSDE\" ' I 




128 


3.023(0.348) 


5.171(0.470) 


512 


2.883(0.176) 


4.973(0.229) 


2048 


2.858(0.094) 


4.936(0.113) 


8192 


2.859(0.053) 


4.938(0.07) 


32768 


2.867(0.025) 


4.948(0.033) 



When t = 



M 


-U,:r,iV,AJ jvj^^ 
^ Q,BSDE\" ) 


Fo'^'^^'^"(cr^'*^) 


128 


3.233(0.349) 


1.884(0.343) 


512 


3.008(0.171) 


1.647(0.161) 


2048 


2.956(0.09) 


1.589(0.078) 


8192 


2.948(0.051) 


1.577(0.047) 


32768 


2.954(0.024) 


1.582(0.023) 



For g2{x,y) = logx + O.by 
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M 


^ tio,BSDE[" ) 




128 


2.925(0.35) 


3.828(0.369) 


512 


2.814(0.178) 


3.706(0.190) 


2048 


2.795(0.094) 


3.685(0.101) 


8192 


2.8(0.054) 


3.690(0.058) 


32768 


2.81(0.025) 


3.701(0.027) 




M 


^ ti5,BSDE\" I 




128 


3.023(0.348) 


5.168(0.403) 


512 


2.883(0.176) 


4.987(0.208) 


2048 


2.858(0.094) 


4.950(0.111) 


8192 


2.859(0.053) 


4.950(0.062) 


32768 


2.867(0.025) 


4.960(0.029) 



When t = 



M 


^ o.bsdeK" ) 


F:^'"'^"(a^^^O 


128 


3.233(0.349) 


1.882(0.309) 


512 


3.008(0.171) 


1.659(0.154) 


2048 


2.956(0.09) 


1.604(0.081) 


8192 


2.948(0.051) 


1.595(0.046) 


32768 


2.954(0.024) 


1.600(0.022) 



For gsiy, z) = O.lz + 0.5y: 
When t is close to maturity 



When t = 



M 


^ ti3,BSDE\" ) 


^ tio y'^ ) 


128 


2.925(0.35) 


3.163(0.403) 


512 


2.814(0.178) 


3.027(0.206) 


2048 


2.795(0.094) 


2.999(0.105) 


8192 


2.8(0.054) 


3.006(0.059) 


32768 


2.81(0.025) 


3.015(0.028) 




M 


^ ti5,BSDE[" I 


F--'-(a-^-) 


128 


3.023(0.348) 


3.541(0.465) 


512 


2.883(0.176) 


3.366(0.231) 


2048 


2.858(0.094) 


3.334(0.113) 


8192 


2.859(0.053) 


3.339(0.069) 


32768 


2.867(0.025) 


3.350(0.033) 




M 


-U,:r,iV,iU , N,M\ 

^ o.bsdeK" ) 




128 


3.233(0.349) 


2.867(0.343) 


512 


3.008(0.171) 


2.654(0.164) 


2048 


2.956(0.09) 


2.603(0.077) 


8192 


2.948(0.051) 


2.595(0.045) 


32768 


2.954(0.024) 


2.6(0.023) 
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We note that we obtain for i = 0, an approximation of the solution of the BSDE which is very 
close to the approximations and the reference price (2.95) in [10]. We can observe, for the same 
parameters, the impact of of the function gi on YiY^^^' ' = 1.58). 

We can also examine the impact of the variable z in the function g on the BDSDE's solution: 
Examining the results in the tables for 171 which depends on z and for g2 which docs not, we 
observe that when g doesn't depend on z, the standard deviation of our solution is smaller, either 
we are close to maturity or at t^. However, for the BDSDE's solution value, when we are close to 
maturity (at iis), it is closer to the BSDE's one when there is dependence of the function g on z. 
At to, the opposite becomes true. 

Finally, we see on the following figures the impact of the function g on the solution; we variate TV, 
M and d as in [11], by taking these quantities as follows: First we fix di = 40 and ^2 = 180. Let 
j e N, we take am = 3, /3 = 1, TV = 2(^\/2)^^-^\ M = 2(a/2)"«(^-i) and d = / {^/2)^'-^^'^P+^^/^ . 
Then, we draw the map of each solution at t = with respect to j . 





Reference price 




BSDE 


— *- 


BDSDEforgI 




BDSDE for g3 




Figure 1. Comparison of the BSDE's solution and the BDSDE's one: The solution of the BSDE is with circle 
markers (in red), the solution of the BDSDE for gi(x, y, z) = O.lz + 0.5j/ + logix) is with star markers (in green) 
and the one for 32(3/, ^) = O.lz + 0.5y is with cross markers (in blue). The reference price for the BSDE is with x 
markers (in red). Confidence intervals are with dotted lines. 



In the second figure, we are interested in analyzing the dependence of the BDSDE's solution on 
the variable z in the function g. So, we variate the parameters iV, M and 5 as above and we draw 
the maps of BDSDE's solution at t = with respect to j for different values of the contraction 
constant a. 



We see clearly from figure 2 that the larger a is, the stronger the impact of z on the solution 
becomes. We finally notice that when a > 0.7 and j > 8, the solution explodes and takes very big 
values. 
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Figure 2. The BDSDE's solution for different values of a which is the contraction constant in z of the coefficient 
9- 



7. Appendix 

7.1. Proof of Proposition 4.2. 

To simplify the notations, we restrict ourselves to the case k = d = I = 1. {DgY, DgZ) is well 
defined and from inequalities (2.4) and (4.1), we deduce that for each 9 <T 



E[ sup \DeYsf]+E[[ \DeZ,\^ds]<C{l + \x\^). 

t<s<T Jt 



(7.1) 



We define recursively the sequence (1"™, Z™) as follows. First we set {Y'^, Z'^) = (0, 0). Then, given 
(y"-i,Z"-i), we define (r",Z'") as the unique solution in S^([t,T]) x H^^^([i,r]) of 

i^" = $(X^")+ / f{r,Xl''\Y^-^,Z'^-^)dr+ I .g(r, X*'", F;""!, Z,""!)^ - / 

J s J s J s 

We recursively show that (y™,Z") G ^^([i, T], pi-^). Suppose that (y",Z'") G B^{[t,T],I])^'^) 

and let us show that (y'»+i^ ^'»+i) ^ B^{[t,T],Ii^'^). 

From the induction assumption, we have $(Xt) + /J" /(r, S™)dr G D^'^. 

We have g{r, S^!') G D^-^ for all r G [t, T]. From Lemma 4.2, we have /^^ g{r, 'i:''^)tBr G D^'^. then 

= i?[ci>(x^--) + /(r, + .9(^ ^T)^r\TZ V -^cr] e 

where S™ := (X;'^,y;", Z™). 
Hence 



^ /(r, + 5(^, )S^ - e 
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It follows from Lemma 4.1 that e Ml^^{[t,T],^}'^^'^) and we have DgY^+^ = DgZ^+'^ = 

ioTt<s<e and for 6* < s < T 

DgY;!'"+^ = V$(X^^)i:»e^^'' (7.2) 

+ / (v,/(r,S;")7?eX, + Vj,/(r,I]™)i^,r™+V,/(r,I]™)-De^r)c?^ 



From inequality (2.4). we deduce that for each 6 <T 

E[ sup l-Derr+T] + E[ r \DgZ^'+'\^ds] < 0(1 + \x\^). 

t<s<T Jt 

It is known that inequality (2.4) holds for Z'"+^) and so we deduce that 

||>^"+'||l,2 + ||^™+'||l.2 <00, 

which shows that (y™+i^^™+i) £ S^([t, T], D^'^). Using the contraction mapping argument as in 
El Karoui Peng and Quenez [9], we deduce that (y™+\z™+i) converges to (F, Z) in S2([f,r]) x 
n^{[t,T]). We win show that {DgY'^.DgZ"') converges to {Y^.Z'^) in L'^{Vl x [t,T] x [t,T],dP® 
dt ® dt), where F/ = Zf = for ah t < s < 6* and (F/, ^'f , 6* < s < T) is the solution of the 
BDSDE. 

r/ = V$(X^"^)i:»e^^'' (7.3) 
+ ^ (v./(r, E,)i?eX, + V,/(r, E,)y/ + V,/(r, 'Zr)zf)dr 

+ j (V.5(^ E,)i^e^, + Vj,5(r, E,)y,'' + V,g(r, E,)Z,^) dS, 

ZldWr. 

From equations (7.2) and (7.3), we have 

DgY:'+^ - r/ = ((V,/(r, E™) - V,/(r, E,))i^eX;'" 

+V^/(r, E::')i^e>;™ - V J(r, E,)K; + V,/(r, E™)i?(,^™ - V,/(r, E,)Z,»)dr 

+ j ((V.g(r, E™) - V,.g(r, E,))i?e^*^- + Vy5(r, E;")i?e>;" - V^.g(r, E,)y,'')dB, 

+ (V.g(r, E::')i^e^™ - V,5(r, E,)Z,«)^ 

- / (7?e^r+i - Z^dWr. 

J s 
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From Proposition 4.1, we have 

E[ sup \DgY:^+^ - y/|2] +E[I \DgZ^^+^ ~ Zl 
e<s<T Js 



dr] 



(7.4) 



< CE 



dr 



+CE 



Therefore, we obtain 



dr 



7t Jt Jt Jt 

< CE[j^ ^ \S^:,\'drd9]+CE[l^ ^ Ipl^el^drdO], 



where 



and 



= (V./(r, E™) - V,/(r, + V,/(r, I]™)^,'' - V,/(r, 

+ V,/(r,I]™)Z,,«-V,/(r,I],,)Z^ 



P":^ = (Vx.9(r, E™) - V,g(r, E,))i^e^;'^ + Vj,5(r, E™)y,^ - Vj,g(r, E,)y/ 
+ V,.9(r, E™)Z,« - V,g(r, E,)Z^ 



(7.5) 



(7.6) 



(7.7) 



From the definition of (,5;';'e)t<^,e<T, we have i;[/f/f|J™g|2drd6l] < C {A,n{0,t,T)+B„,{0,t,T))de, 
where 



A„(0,<,T) 



E 



|(V,/(r, S-) - V,/(r, i:,.))DeXl'^\'dr 



Bmie,t,T) = E 1^ |(V,J(r,E,)-V,/(r,E™))y/pdr 



E 



(V,/(r,S,)-V,/(r, E™))Z« rdr 



Moreover, since Vxf is bounded and continuous with respect to (x, y, z), it follows by the dominated 
convergence theorem and inequality (2.3) that 



lim / A^{e,t,T)d0 = 0. 



(7.8) 



Furthermore, since Vj,/ and Vzf are bounded and continuous with respect to {x,y,z), it follows, 
also, by the dominated convergence theorem and inequality (2.4) that 



lim / B^{e,t,T)de ^0. 

m—>-oo J ^ 



(7.9) 
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From the definition of (p™(,).<r,e<T, we liavc E\^^ |p™epdrd6'] < C (A'„(6', t, T)+B',^{e, t, T))d0 
witli 



E 



+ E 



Similarly as shown above, since V^g and V zg are bounded and continuous with respect to (x, y, z) 
we can show that: 



lim / A'^{0,t,T)de = lim / B',^{e,t,T)d0 ^ 0. 
Plugging (7.8), (7.9) and (7.10) into inequality (7.5), we deduce that 

lim E[[ I \DgY^'"+^ -Y,^\^dsde]+E[[ [ \DgZl''+'^ ~ Z^J^dsdO] = 0. 



(7.10) 



It follows that (F™, Z") converges to {Y, Z) in L^{[t, T], 
is given by (Y^ , Z^) which is the desired result. 



a, 2 ^ inii,2 



) and a version of {DgY, DgZ) 



□ 



7.2. Second order Malliavin derivative of the solution of BDSDE's 

We apply know similar computation to get the second order Malliavin derivatives representation 
of the solution of BDSDE 's, so we will omit the proof. 

Proposition 7.1 Assume that (H2)-(H3(ii))-(H4(ii)) hold. We fix t e [0,T]. Then for each 
t<0 <T, {DeY,DgZ) belongs to ^^([t, T], D^'^). For each t < v < T and 1 < i, j < d, 

DiD'gYs = DiDlZl' = 0, 1 < n < d, ifs<9\/v, 

and a version of {DlD^gYs, DlDgZs)v\je<s<T is the unique solution of the equation: 

DiDlY, = Ti($) + T2(/) + n{g) + mw), 

where 

k 

ni — 1 
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T2U) = r ^ (v,((v,/^H^^*'^>^r.,^.))^ix;'"(^^^;■"r 

ni = l 



^3(5) 



V,f{r,Xy\Yr,Zr)DlDlXl^-)dr 



I ( E ^y{{^yfr{r,Xr.Yr,Zr))DlYr{DlYrf 



n\—l 



Vyf{r,X'^--,Yr,Zr)DiDiYr)dr 

d A; 



12 = 1 ni=l 
d „T 

1 s 



"2 = 1 



J2 r ^Ai^x9"'rHr,Xl^\Yr,Zr))DiX*'^DlXl^-r'^' 

n3 = l-^*' ni = l 
ri3 = l"^* ni = l 

E E r E ^--((^--5"^r(^'^*'^^-^'0M^r(^e^rrSs^ 

n3^1 712—1 ^ Tl-i— 1 

E E f^^-^9'''ir,X','^,Yr,Zr)DiDlZ;^^^^, 



"3 = 1 "2 = 1 



T4(iy) = - E f DiDlZl'-dW:^\ 

"2 = 1 



(-^"')i<i<'i denotes the j-th column of the matrix z, {g'^^)i<na<i denotes the n^-th column of the 
matrix g, B = {B^ , B''), (D^X*>=^)"i is the ni-th component of the vector (D^X*'=^), (f^Fr)"' 
is the ni-th component of the vector (DgYr) and (DgZ!^^)""^ is the ni-th component of the vector 

7.3. Some estimates on the solution of the FBDSDE 

Lemma 7.1 Let {b^,a^) and {b'^,a'^) he the standard parameters of the SDE (2.1) with initial 
condition x^ (resp. x^). We assume that (HI) holds. Put 6Xs = X^ — X^, Sbg = {b^ — b'^){Xl) 
and 5a, = {a^ - a^){Xl).Then 

\\X%.<C{l + \xf). 
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For all Si,.S2 S [0,T], we have 
E 



sup \Xl^Xl^\ <C(l + |xn|s2-si|, 

Si<U<S2 



and for all si < s < S2, we have 



\\sx\\ 



< c 



„2|2 



\S2-S,\+E[ f ' \dbs\^ + \d<Js\''ds]), 
J si 



where C is a generic constant depending only on K, T, (6^(0), cr^(O)) and (6^(0), cr^(O)). 

Lemma 7.2 Let {X, Y, Z) be the solution of the FBDSDE (2.1)-(2.2). We assume that (H1)-(H2) 
hold. Then, we have 



|y|u + ||z||n.„ <c{i + \x\^). 



id also 



E 



sup 



iFn-l^tp] <c[{l + \x\'')\s-t\+ I iZ.pdr) 



(7.11) 



(7.12) 



Proof. The technics used to proof these estimates are classical in the BSDE's theory so we omit 
it. 
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